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The Sivers distributions recently extracted from semi-inclusive deep inelastic scattering data [lj 
are used to compute estimates for Sivers asymmetries in Drell-Yan processes which are being planned 
at several facilities (RHIC, COMPASS, J-PARC, PAX, PANDA, NICA (JINR) and SPASCHARM 
(IHEP)). Most of these asymmetries turn out to be large and could allow a clear test of the predicted 
sign change of the Sivers distributions when active in SIDIS and Drell-Yan processes. This is 
regarded as a fundamental test of our understanding, within QCD and the factorization scheme, of 
single spin asymmetries. 
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. I. INTRODUCTION 

' The experimental study and theoretical understanding of transverse Single Spin Asymmetries (SSA) has been, and 
still is, one of the most challenging issues in high energy hadron physics. The original, widespread opinion that 
O ■ hadronic SSAs should be negligible in any high energy process, due to the simple, helicity conserving pQCD and 
CD | Standard Model elementary dynamics 0], has been proven wrong in a great number of cases, actually in most of 
i-C ■ the spin measurements so far performed. It often happens that the same approach, based on the collincar QCD 
factorization scheme at leading twist, which successfully describes unpolarized data, cannot explain the large spin 
effects observed in the same kinematical regions. 
>• The attempts of explaining the data and reconciling them with pQCD dynamics have produced a much deeper 

understanding of the nucleon structure, involving partonic intrinsic motion, and much progress in clarifying the 
mechanism, at the partonic level, responsible for providing the phases and helicity flips necessary for a transverse 
SSA. Many issues, like the universality of these mechanisms and their correct insertion into a factorized scheme, are 
still under debate and need further investigation. 

We consider here the so-called Sivers asymmetry [3, 4| , related to the intrinsic motion of partons inside a transversely 
polarized proton, according to the distribution 
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f q/p i(x,k±) = f q/p (x,k±) + - A N f q/p r{x,ki_) S- (P x k_i_) (1) 

= f q/p (x, k ± ) - ^ k±) S-(Pxk ± ), (2) 

' nip 

which gives the number density of unpolarized quarks q (or gluons) with intrinsic transverse momentum k± inside 
a transversely polarized proton p\ with three-momentum P and spin polarization vector S; A f q / p i(x,k_i), or 
f^(x,k±) following a different common notation [j|, is the Sivers function. A coupling of the transverse motion 
of unpolarized quarks and gluons to the nucleon spin can only be related to their orbital motion, and one expects 
that the observation of the Sivers asymmetry signals partonic orbital angular momentum 0, 0]. Thus, the Sivers 
distribution is of very special interest. 

A clear observation of a non zero Sivers distribution has been obtained by the HERMES collaboration in Semi- 
Inclusive Deep Inelastic Scattering processes (SIDIS) @,[j|. A similar measurement by the COMPASS collaboration 
has given a result compatible with zero EL EH ; however, such a result was obtained on a polarized deuteron target 
— rather than a hydrogen one, as for HERMES — and the small result could be explained as a cancellation between 
opposite contributions from u and d quarks. That has allowed the extraction, from data, of the Sivers distribution 
function [l . Il3l . fli , ITU Il6j . Large SSAs observed in p^p — ► n, K + X inclusive processes at Fermilab [13, EM HH and 
RHIC [2(1 2l| could also be explained as a manifestation of the Sivers effect [22|, HH, [24| . A recent preliminary result 
from COMPASS operating on a transversely polarized proton target [2o| and still compatible with zero, might, if 
confirmed, be of difficult interpretation. 
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From the theoretical point of view the Sivers distribution function 0, [|| had a difficult early life. For some time 
its very existence was much debated and, despite its phenomenological success [22|,[23|, fundamental parity and time- 
reversal properties of QCD seemed to forbid it j2^|. Some attempts were made of invoking initial state interactions [23| 
or unusual time-reversal properties 27 1. An explicit model calculation of the Sivers effect in SIDIS [28| much clarified 
the situation, showing the crucial role of final state interactions. The analysis of Ref. [2(| was then reconsidered [29[ 
and the proof that the Sivers asymmetry should vanish because of time-reversal invariance was shown to be invalidated 
by the path-ordered exponential of the gluon field in the gauge invariant operator definition of parton densities. Thus, 
the gauge links (gluon exchange in final state interactions) seem to allow a consistent QCD picture of the Sivers effect. 
The same analysis [29[ which led to the conclusion that the Sivers distribution has a full right to existence in QCD, 
also showed that time-reversal properties imply that the Sivers asymmetry must be reversed in sign when acting in 
SIDIS and Droll- Yan (DY) processes. This was explicitly confirmed by the extension of the model for Sivers effect in 
SIDIS to the DreU-Yan case [H. 

The experimental measurement of the Sivers SSAs in SIDIS and Drell-Yan processes, and the observation of the 
sign change of the Sivers distribution, is one of the most important tests of our understanding of the origin of SSAs in 
QCD and the factorization scheme. In this paper we exploit the Sivers functions recently extracted from SIDIS data 
change their signs, and give estimates for Sivers SSAs in Drell-Yan processes; these asymmetries can be measured 
in several experiments, to be performed at existing facilities (COMPASS, RHIC), developing (J-PARC) or proposed 
ones (PAX, PANDA, JINR, IHEP). 

An early discussion of the contribution of the Sivers function to SSAs in Drell-Yan processes, for RHIC experiments, 
was presented in Ref. (3l| ; information on the Sivers distributions was obtained from data on SSAs in inclusive hadronic 
production, nj_p — > it X, rather then from SIDIS data. This paper updates and extends subsequent studies performed 
in Refs. [3, [HI and [33|. In Ref. (3] predictions for Sivers SSAs in Drell-Yan processes at RHIC and PAX, based on 
a first extraction of the Sivers functions from SIDIS data, with no sea contributions, were given. A similar procedure, 
for RHIC experiments, was followed in Ref. [l5[. In Ref. [32| estimates for SSAs at RHIC, PAX (in fixed target mode) 
and COMPASS were given based on a simple expression of the (opposite) u and d Sivers functions, extracted from a 
two-parameter fit of the SIDIS HERMES data, and combined with some models for the antiquark Sivers distributions. 
Here, we use the recent [l|, much more structured, valence and sea, Sivers functions obtained from a fit of SIDIS 
asymmetries for pion and kaon production, measured by HERMES and COMPASS collaborations. We give estimates 
for the full set of experiments which cither have been proposed or are being discussed in many laboratories worldwide. 



II. FORMALISM FOR THE SIVERS EFFECT IN DRELL-YAN PROCESSES 

We consider the leading order (LO) parton model cross section for Drell-Yan processes, B — > £ + £~X (A* = 
; B — p,p, 7r), in the hadronic cm. frame, in which one only observes the four-momentum of the final £ + £~ pair, or 
related quantities: 

q = (<Zo,<?T>9i) q =M y =- In x F = —= s = {p A +PB) • (3) 

2 qo - qL V s 

In the kincmatical region 

q% < M 2 k ± ~q T , (4) 

where qt = \Qt\ an( ^ ^J- = I^jJ is the magnitude of the parton intrinsic motion, at 0(k±/M) the Sivers SSA is simply 
given by [3l| 

foAiB^e+t-x _ da Ai B -.i+e-x _ dgT _ da l 

da AiB^i+e-x + da AiB^i+e-x ~ d(J i + dcr i ( 5 ) 
J2 q J d 2 k ±1 d 2 k ±2 S 2 (k ±1 + kj_ 2 - q T ) A N f q/Ar (xi,kj_i) fq/ B {x2,k 1 _ 2 ) 



2 T iq I d 2 k_n d 2 k ±2 S 2 (k ± i + fc ±2 - q T ) f q / A (xx,k±i) fq/ B {%2,k± 2 ) a™ 



(6) 



with 



-_ W M ±v _ ±x F + ^/xl+4Mys 

9AP 1,2 ~ J~s 2 [n 



where the subscripts 1 and 2 label the quark and antiquark participating in the elementary LO scattering process, 
qq — v £ + £~ . Notice that, with the definition of xf adopted in Eq. ([3]), one has 

M 2 

xf = x\ — X2 \xf\ < 1 • (8) 
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The sum in Eq. (|6|) runs over all quarks and antiquaries {q = u, u, d, d, s, s) and da stands for 

d A o~ 1 d A o~ d A a 1 d A a 

dydM 2 d 2 q T s dxidx 2 d 2 q T ' dx f dM 2 d 2 q T 2 d*q 

Notice that in obtaining Eq. §§§ one uses (see Eq. (Q])): 

f q /AT(x,k±) + f q/A i(x,k ± ) = 2f q/A (x,k±) (10) 
fq/AT (x, k±) - f q/A i (x, k±) = A N f q/Ar (x, k±) S-(Pxk ± ) = A N f q/Ar (x, k±) . (11) 

The SSA © depends directly on the Sivers functions A N f q / A i . 

We use here the same factorized expression of the cross sections which holds in the collinear configurati on, g eneraliz- 
ing it to the case of unintegrated, or transverse momentum dependent (TMD), partonic distributions [24l. l33l| A most 
general treatment of unpolarized and polarized Drell-Yan processes in such a scheme has very recently appeared }34| . 
Factorization for SIDIS and Drell-Yan processes in the kinematical regime we are considering here, Eq. has been 
proven in QCD (3f| [H, H?]], resulting in the same parton model TMD-factorization adopted here and in Ref. [34[, with 
the addition of an extra soft factor S 1 , which takes into account transverse motion originating from soft gluon emission 
(see, for example, Eqs. (40) and (41) of Ref. [H[). Such a factor gives an (unknown) additional contribution, both to 
the numerator and denominator of Eq. (JSJ), of order a s , and is neglected here. The TMD partonic distributions we 
use are obtained by fitting experimental data and take into account all sources of intrinsic motion. 



A. About sign and azimuthal angle conventions 

As the issue of the sign of the Sivers asymmetry in Drell-Yan processes is so important, let us discuss in details the 
choices adopted here and their motivation. We define our kinematical configuration with hadron moving along the 
positive z-axis, and hadron B opposite to it, in the A — B center of mass frame. We choose the "up" (|) polarization 
direction as the positive y-axis (0s = t/2). The transverse momenta have azimuthal angles 

q T = 9t(cos0 7 , sin0 7 , 0) k ±i = fcj_ 4 (sin ^ , cos<^, 0) (z = 1,2), (12) 

so that the mixed product S ■ (P x fej_i) gives an azimuthal dependence sin(</>s — ifi) = cos tpi which, upon integration 
on fej_x, yields a sin(0s — 7 ) = cos</> 7 dependence of the Sivers asymmetry (see Eq. (|^D|) below). 

Notice that, contrary to the usual study of angular dependences of Drell-Yan processes which is mainly performed 
in the so-called Collins-Soper reference frame [39(, we work in the hadronic cm. frame. At least for the purpose of 
studying the Sivers asymmetry, this frame is much more convenient and directly related to experimental measurements. 

In order to collect data at all azimuthal angles, following what is usually done in semi-inclusive deep inelastic 
scattering processes, both experimental results and theoretical estimates can be discussed for the azimuthal moments 
of the asymmetry. We follow Refs. [13, E3] and [l5| and choose as a weight the sin(0 7 — <ps) phase. We then have: 

A sm(0 T -0 s ) = Jo* #7 [ rffjT ~ d(ji \ Sin (^7 ~ <h) ^ 

N 3 Jo* <#7 + dai] 



J d</> 7 


'E q Jd 2 k^d 2 k ±2 S 2 (k ±1 +k^ 2 


- Qt) & N fq/Ai (xi,kxi) fq/ B {x2, k ±2 ) a% q 


sin 


(07 ~ ^S) 


/ d</> 7 


'J2 q jd 2 k^d 2 k ±2 6 2 (k^ 


_i +k ±2 - q T ) f q /A(xi,kxi) fq/B{x 2 ,kj_ 2 ) 







Notice that such a choice, combined with the sin(0g — 7 ) dependence associated with the Sivers function [see 
comments following Eq. (Ti"2|) ] implies an overall [— sin 2 (0 7 — <fis)] factor in the numerator of Eq. (fl3|) . 

The above asymmetry is, in general, a function of xf (or y), M and qT- In the sequel we shall study it as a function 
of one variable only, either xf or M; we will always integrate both the numerator and denominator of Eq. (|13[) over 
qr — covering the range in which the factorized approach with unintegrated distribution functions is supposed to 
hold (as detailed in the captions of Figs. [5]© — and on one of the remaining variables according to the kinematical 
ranges of the corresponding experiments. 

The case in which the polarized hadron moves along — z, that is the process B A^ — > £ + £~X, deserves a special 
comment. In such a case — keeping the same definition of f polarization and the same reference frame — the Sivers 
mixed product has an opposite sign (due to the opposite A^ momentum) and yields a Sivers asymmetry proportional 
to sin(0 7 — 0g). In this case the overall factor in the numerator of Eq. (TIB")) is [+sin 2 (0 7 — 0s)]. This agrees with 
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what is done in SIDIS processes, £p^ — > £hX, in the 7* — p^ cm. frame Q. To summarize, we shall give estimates 
for the quantities: 

A f^-*s) {A i B ^ rx . a; Fj M, qr) = -Af^-'-^iBA^^X- -x F ,M,q T ) . (14) 
The equality holds due to rotational invariance. 



III. ESTIMATES FOR FORTHCOMING EXPERIMENTS 



In order to give estimates for the Sivers asymmetries in Drcll-Yan processes — and test the crucially important 
sign change when going from SIDIS to DY — we only need to insert the Sivers functions extracted from the analysis 
of SIDIS data into Eq. (fl~3|) . We use the results obtained in Ref. [H, which adopted a Gaussian factorized form both 
for the unpolarized distribution functions: 



,-fci/<fci> 



fg/ p (x,k±) = f q (x) 

and for the Sivers distributions: 

A N f q/p r(x,k±) = 2Af q (x)h(k±)f q/p (x,k±) 

= A N f q/p i(x)h(k ± ) 

where 



(ki) = 0.25 GeV 2 



(15) 



(16) 



M q (x)=N q x*«(i-xf« {aq+ :± jq — 

h(k±) 



(17) 
(18) 



The values of the 11 best fit parameters N q (q = u, d, s, u, d, s), a q (q = u, d, sea), (3 (same for all q) and M\ can be 
found in Table I of Ref. [l|, where their uncertainty is also explained in details. 

Notice that the above factorized expressions allow, at 0(k±/M), an analytical integration of the numerator and 
denominator of Eq. (fill)) , resulting in 



.4 



i(<t>-r—<t>s) 



N 



(x F ,M,q T ) 



J d<p 7 [N(x F , M, qr, 7 )] sin(</> 7 - </>g) 
/ d^ [D(x F ,M,q T )] 



(19) 



with (see Eq. ©): 

N(x F ,M,q T ,<t> 7 ) 
Ana 2 



d 4 a l 



dx f dM 2 d 2 q T dx F dM 2 d 2 q T 



E 



q T (k 2 ) 2 ex P [-q 2 T /((k 2 ) + (k 



2 ) 

±21 



sin(0 s - 7 ) 

(20) 



and 



D(x F ,M,q T ) 



d 4 al 



d 4 a l 



dx F dM 2 d 2 q T dx F dM 2 d 2 q T 



d 4 a unp 



in a 2 
9M 2 s - 



E 



fq/A{xi) fq/B(x 2 ) 



dx F dM 2 d 2 q T 
exp[-^/((fcj 1 ) + (fci 2 )' 

7T[(fci 1 )+(fci 2 )] 



(21) 



Notice that we have defined 



M 2 



(22) 
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and that xi,X2 have the values given in Eq. (J7J). Although in the above equations we have allowed for the possibility 
of having different values of (fcj^ 2 ) f° r hadrons A and B, in our numerical evaluations we shall take them equal, as 
given by Eq. (|15j) , even when one of the hadrons is a pion. 

Our estimates are obtained using, for the proton, the partonic distribution functions (PDF) of Ref. [4l[ (GRV98LO, 
consistently with the choice adopted in Ref. [l[) and, for the pion, those of Ref. [42|; all the PDFs are evolved (at LO) 
at the Q 2 = M 2 scale of interest. We have checked that a different choice for the PDFs of the pion [43| yields very 
similar results for the 7r~ and results in agreement with those shown here, within the uncertainty bands, for the 7r + . 
Our results are presented in Figs.QJH] The shaded bands reflect the statistical uncertainty in the extracted values of 
the Sivers functions, as explained in Ref. [![. Let us comment on each figure. 

• Fig. [TJ Sivers functions. 

For completeness, the Sivers functions used in this analysis are presented. The sign has been reversed with 
respect to the Sivers functions obtained from SIDIS data in Ref. according to the findings of Refs. (29l. [30j: 

AV = -A N f sims . (23) 

The first moment of the Sivers function, left panel, is defined as: 

A N f^ (*) ^ Jd 2 k^^- A N f q/p1 (x, fe x ) = -f^ 9 (x) . (24) 

This quantity is, up to a constant resulting from the d 2 k± integration, the same as the function A f q /Ai (x±) 
appearing in Eqs. (Tj"6|) . ([20)) . 

• Fig. H COMPASS. 

as a function of xf = x± — X2, left plot, and of M, central plot, for the COMPASS planned 
experiment, ir^ p^ — » n + /i~X. We have integrated both the numerator and denominator of Eq. (|19p over qp 
in the range (0 < qp < 1) GcV, which is within the region of validity of our approach. The other integration 
ranges are (4 < M < 9) GeV, at fixed xf, for the left plot and 0.2 < xf < 0.5, at fixed M, for the central plot. 
The pion beam energy, in the laboratory frame, is taken to be 160 GeV, corresponding to = 17.4 GeV. 

In the right panel we show the kinematical region covered by the COMPASS experiment, i.e. the range of allowed 
X2 values, which refer to the polarized proton distributions, as a function of xf- It is interesting to notice that 
positive values of xf correspond to (or at least overlap with) the x region explored by the SIDIS experiments 
(x < 0.35) and the data used to extract our Sivers functions. Instead, negative values of xf correspond to 
larger X2 values, a region yet unexplored by other experiments, so that our estimates arc based on the assumed 
functional form of the Sivers function — not constrained, in that region, by any SIDIS data — and are bound 
to have much larger uncertainties. 

Finally, it is important to remark that, as the COMPASS experiment will involve charged pion beams, ir~ (ud) 
and Tr + (ud), the dominant elementary process contributing to the asymmetry will be u^-u p — > p. + fi~ for the 
7T _ beam and d n +d p — > for the 7r + beam. Consequently, the prediction of the sign change of the Sivers 

function in SIDIS and Drell-Yan processes [1^, [3(| can be clearly tested, as the sign of the asymmetry for tt~ 
(tt + ) beam is given by the sign of the u (d) Sivers function, which is well established. 

• Fig.d COMPASS (low M). 

We show our estimates for the single spin asymmetry and the corresponding allowed kinematical region in the 
X2 — xp plane, for COMPASS experiments performed at lower M values, (2.0 < M < 2.5) GeV and in a more 
limited qr region, (0 < qr < 0.4) GeV (so that the constraints of Eq. (0} remain true). This exploits the 
possibility of DY measurements in the low mass region, below the J/ip and above the p resonance peaks, where 
data should be much more abundant. 

• Fig. |H RHIC. 

^sm(0 T <t> s ) j g s i 10wn as a function of xp and M for RHIC experiments, p^p — » £ + £~X, at y^s = 200 GeV. The 
integration ranges for qp and M are the same as in Fig. [21 with the further constraint < y < 3, according to 
the experimental kinematical cuts. The right panel shows the kinematical region of RHIC for x\ as a function of 
xp. This region covers the range already explored by SIDIS measurements (x < 0.35) where the Sivers functions 
are reliably constrained, and expands into much higher values of x. The maximum value of the asymmetry (~ 
10%) shown in the left panel is expected at xf — 0.2 which corresponds (see the right panel) to x\ ~ 0.2, where 
the valence Sivers functions reach their maximum, see Fig. [1] 
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The asymmetry for larger values of xf ( > 0.4) is obtained by using our extracted Sivers functions at large values 
of x where SIDIS data offer no constraint. This reflects into the huge uncertainty band. Let us remind that in 
Ref. [I[ we have assumed the (3 parameters, which determine the large x behaviour of the Sivers functions, to be 
equal for all flavours. With such an assumption the asymmetry decreases fast for high values of xf, although 



with a large uncertainty. The parameterization of our previous analysis 14j , where (3 U and (3d were different for 
u and d flavours, would lead to a larger asymmetry at high xf (see Fig. 6 of Ref. Measurements in the 

region of high xp could then offer an opportunity to test the flavour-blind (3 assumption. Moreover, data in the 
negative xf region would test the contribution of the sea Sivers functions, as first pointed out in Ref. [32| . 

• Fig. U PAX. 

Estimates are given for the planned GSI-PAX experiment [Hf, p^p — > £ + £~X, in the asymmetric collider mode 
at ^/s = 14.14 GeV, with a polarized proton beam. This experiment would have the big advantage of having 
proton and antiproton beams, thus increasing the number of Drell-Yan events. is shown as a function 

of Xf (integrating over (4 < M < 6) GeV, with the constraint \y\ < 1), and as a function of M (integrating over 
xp with \y\ < 1); in both cases qp is integrated in the range (0 < qp < 1) GeV. 

The negative xf region gets most contribution from the valence x\ region, which is well explored by SIDIS 
experiments and where the Sivers distributions are reliably known. In this region the asymmetry is large (10- 
15% in size) and sensitive only to the u and d Sivers functions. On the contrary, at large and positive xf values 
the situation becomes similar to that described for RHIC, as we need the Sivers functions in the large x-range, 
which are more poorly known. Once again, this results in large uncertainties of our estimates. 

The PAX collaboration aims at having a beam of polarized antiprotons, which, in combination with a beam of 
polarized protons, would allow a unique direct way of measuring the transversity distributions In such a 
case, one could also measure the single spin Sivers asymmetry in the process pp* — > £ + £~X; this is related to 
the one considered here by charge conjugation invariancc and the relation of Eq. (|14|) . 



A s ^'^ s \pf - 7 **; x F ,M,q T ) = -Af^ s \plp~>rX; -x F ,M,q T ) . (25) 

• Fig. H PANDA. 

The GSI- PANDA experiment [45[ will run with an antiproton beam scattering off a proton target, at the 
energy Ep = 15 GeV, corresponding to yfs = 5.47 GeV. It might be worth trying to have a polarized target, 

pp^ — ► £ + £~X. We then show predictions for A^ 1 ^ ^ s " > as a function of xf] the integration ranges, due to 
the moderate energy, are (2.0 < M < 2.5) GeV and (0 < qp < 0.4) GeV. The kinematical range which could be 
explored is entirely overlapping with that covered by the SIDIS data, so that this experiment would supply a 
perfect consistency check of the HERMES and COMPASS results, and of the crucial Eq. (|23|) . The asymmetry 
turns out to be sizable and well definite in sign over the full kinematical range. 

• Fig. J-PARC. 

J-PARC might measure Drell-Yan single spin asymmetries A^ 1 ^" 1 generated in polarized proton-proton 
collisions at E p = 50 GeV (HI, corresponding to ^fs = 9.78 GeV (see the figure caption for further information). 
As shown on the right panel, the kinematical range covered by J-PARC would be almost entirely complementary 
to that explored in the SIDIS experiments. These measurements could then be of great importance to obtain 
information on the large- a; behaviour of the Sivers distribution functions. 

• Fig. H NICA. 

NIC A is a planned Drell-Yan experiment to be performed at JINR in Dubna [47|- We present our estimates 
for an unpolarized proton beam scattering off polarized protons at y/s = 20 GeV. As shown in the right panel, 
the region of large and negative xp ( < —0.5) corresponds to the range of x^ values for which we have no 

information from SIDIS experiments. Therefore the measur ements of Af^'^ in this region would provide 
information on the large x behaviour of the Sivers functions. Our results have smaller uncertainties in the 
region —0.4 < xf < 0.1 which corresponds to the X2 valence region and where the asymmetry is related to the 
convolution 4 f a <£> A N f u + f 3 ® A N f d . The comparison between the left and the central plots shows that the 
positive xp region is dominated by the Sivers sea contribution. More information can be found in the figure 
caption. 



Fig. M SPASCHARM. 

The SPASCHARM experiment [48| at IHEP in Protvino (Russia) plans to measure Drell-Yan processes both 
in polarized proton-proton scattering at PLab = 60 GeV, corresponding to ^/s = 10.7 GeV (upper panels) and 
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in pion- proton scattering at E v = 34 GeV, corresponding to yfs — 8 GeV (lower panels). SPASCHARM would 
mainly collect data in the range of low values of M [(2.0 < M < 2.5) GeV] and in the valence region of the 
Sivers functions (see the figure caption for more details). Therefore, SPASCHARM measurements of pion-proton 
asymmetries, like the analogous COMPASS case, could provide a stringent test of the sign change of the Sivers 
function in Drell-Yan and SIDIS processes. 

IV. COMMENTS AND CONCLUSIONS 

We have given clear and simple estimates of Sivers SSAs in Drell-Yan processes, which could be measured in 
experiments which are either in preparation or being planned. Our results show large values of these asymmetries, 
which, despite the uncertainty in magnitude related to the extraction of the Sivers functions from SIDIS data, have an 
unambiguous sign. This sign originates from the current understanding of the origin of SSAs at the partonic level (29| 
which predicts that the Sivers distributions must enter with opposite signs in SIDIS and Drell-Yan processes. Then, 
the mere measurement of this sign would, alone, provide a stringent and important test of the present interpretation 
of the observed SSAs in terms of TMD distribution functions, coupled to initial or final state QCD interactions. 

This test, if successful, would be a significant and decisive step towards a basic and consistent description of 
subtle phenomena like the challenging SSAs in terms of elementary pQCD dynamics and fundamental properties of 
the nucleon structure involving quark intrinsic motion. Drell-Yan processes at small values, exploring the quark 
content of the nucleon without the complications of fragmentation processes, might play a crucial role in the future 
exploration of the proton structure. 
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FIG. 1: The Sivers distribution functions for u, d and s flavours as determined by our simultaneous fit of HERMES and 
COMPASS data in Ref. The sign was reversed according to the prediction of Refs. [2§, 30]. On the left panel, the first 
moment, x A N f^ (x), is shown as a function of x at Q 2 = 2.4 GeV 2 for each flavour. On the right panel, the Sivers distribution, 
x A N f(x,k±), is shown as a function of k± at a fixed value of x = 0.1 for each flavour. In each plot, the highest and lowest 
dashed lines show the positivity limits \A N f\ = 2/. 
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FIG. 2: The single spin asymmetries A"™^' 1 for the Drell-Yan process n ± p l — > fi + fj,~ X at COMPASS, as a function of 
if = X\ — X2 (left panel) and as a function of M (central panel). The integration ranges are (0 < qr < 1) GeV, (4 < M < 9) GeV 
and 0.2 < xf < 0.5. The results are given for a pion beam energy of 160 GeV, corresponding to y/s = 17.4 GeV. The right 
panel shows the allowed region of X2 values as a function of xf- 




FIG. 3: The single spin asymmetries A a ^^~" for the Drell-Yan process n ± p^ — > X at COMPASS, as a function of 

xf = x\ — X2 (left panel). The integration ranges are (0 < qr < 0.4) GeV and (2.0 < M < 2.5) GeV. The results are given for 
a pion beam energy of 160 GeV, corresponding to = 17-4 GeV. The right panel shows the allowed region of X2 values as a 
function of Xf- 




FIG. 4: The single spin asymmetry A 3 ^ 1 ^'' for the Drell-Yan process p^p — > X at RHIC, as a function of xf (left 

panel) and M (central panel). The integration ranges are (0 < qr < 1) GeV and (4 < M < 9) GeV, with the further constraint 
< y < 3. The results are given at y/s = 200 GeV. The right panel shows the allowed region of xi values as a function of xf. 
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FIG. 5: The single spin asymmetry A^" 1 ^ s ' > for the Drell-Yan process p^p — > n + fi~ X at PAX, as a function of xf (left 
panel) and M (central panel). The integration ranges are (0 < qr < 1) GeV and (4 < M < 6) GeV with the further constraint 
\y\ < 1- The results are given at y^s = 14.14 GeV. The right panel shows the allowed region of x\ values as a function of xf- 
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FIG. 6: The single spin asymmetry A a ^^'' for the Drell-Yan process pp^ — > X at PANDA, as a function of xf- The 

integration ranges are (0 < qr < 0.4) GeV and (2 < M < 2.5) GeV. The right panel shows the allowed region of X2 values as a 
function of xf- The results are given at Ep — 15 GeV, corresponding to \fs = 5.47 GeV. 




FIG. 7: The single spin asymmetry A B ^' I '~' ^ for the Drell-Yan process p^p — > X at J-PARC, as a function of xf- The 

integration ranges are (0 < qr < 1) GeV and (4 < M < 5) GeV. The results are given at E p = 50 GeV, corresponding to 
= 9.78 GeV. The right panel shows the allowed region of x\ values as a function of xf- 
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FIG. 8: The single spin asymmetry A"^ for the Drell-Yan process pp 1 ' — > X at JINR (NICA experiment), as a 

function of xf- The integration ranges are (0 < qr < 1) GeV and (4 < M < 9) GeV. The results are given at y/s = 20 GeV; 
on the left panel we include the Sivers sea contributions, while on the central panel we show the results we obtain by including 
only the Sivers valence functions. The right panel shows the allowed region of X2 values as a function of iy. 
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FIG. 9: The single spin asymmetries A a ^ for the Drell-Yan process pp 1 ' — > fi + fi~ X (upper- left panel) and np 1 ' — ► fi + fi~ X 

(lower-left panel) at IHEP (SPASCHARM experiment), as a function of xf- The integration ranges are (0 < qr < 0.4) GeV and 
(2 < M < 2.5) GeV. The energies of the two experiments are, respectively, PLab = 60 GeV, corresponding to y^i = 10.7 GeV, 
and En = 34 GeV, corresponding to y/s = 8 GeV. The right panels show the allowed region of X2 values as a function of xf- 



